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ABSTRACT 

r— i' We present the results of a series of cosmological iV-body simulations of a Vector Dark En- 

\ ergy (VDE) model, performed using a suitably modified version of the publicly available 

. GADGET-2 code. The setups of our simulations were calibrated pursuing a twofold aim: 1) to 

analyze the large scale distribution of massive objects and 2) to determine the properties of 
*^ halo structure in this different framework. We observe that structure formation is enhanced 

in VDE, since the mass function at high redshift is boosted up to a factor of ten with respect 
Q ' to ACDM, possibly alleviating tensions with the observations of massive clusters at high red- 

s'^ , shifts and early reionization epoch. Significant differences can also be found for the value of 

■ the growth factor, that in VDE shows a completely different behaviour, and in the distribution 
& \ of voids, which in this cosmology are on average smaller and less abundant. We further stud- 
ied the structure of dark matter haloes more massive than 5 x 10 13 /i _1 Mq, finding that no 

' substantial difference emerges when comparing spin parameter, shape, triaxiality and profiles 

\ of structures evolved under different cosmological pictures. Nevertheless, minor differences 

i can be found in the concentration-mass relation and the two point correlation function; both 

0^ ' showing different amplitudes and steeper slopes. Using an additional series of simulations 

, of a ACDM scenario with the same Hm and eg used in the VDE cosmology, we have been 

able to establish whether the modifications induced in the new cosmological picture were 
ly-j | due to the particular nature of the dynamical dark energy or a straightforward consequence 

C^i . of the cosmological parameters. On large scales, the dynamical effects of the cosmic vector 

£Nl ' field can be seen in the peculiar evolution of the cluster number density function with red- 

t-H \ shift, in the shape of the mass function, in the distribution of voids and on the characteristic 

■ form of the growth index -f(z). On smaller scales, internal properties of haloes are almost 
• ?-H ' unaffected by the change of cosmology, since no statistical difference can be observed in the 

, characteristics of halo profiles, spin parameters, shapes and triaxialities. Only halo masses and 

5-j ■ concentrations show a substantial increase, which can however be attributed to the change in 

the cosmological parameters. 
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1 INTRODUCTION 



During the last twelve years, a large amou nt of cosmologi- 
cal high precision data on Supernovae la (see iRiess et al I ll998l ; 



IPerlmutter et al.lll999t iGuv et al. 


|2010|) cosmic microwave back- 


ground anisotropies (Larson et al. 


201 ltlSherwin et al.l201 1 ), weak 



lensing jHuterejboiol) , baryon acoustic oscillations (iBeutler et al.l 
1201 lh and large scale structure surveys dAbazaiian et al.ll2009l) has 
provided evidence that the Universe we live in is of a flat geometry 
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and undergoing an accelerated expansion. These observations mo- 
tivate our belief in the existence of an ubiquitous fluid called dark 
energy (DE) that by the exertion of a negative pressure, counters 
and eventually overcomes the gravitational attraction that would 
otherwise dominate the evolution of our Universe. The simplest ex- 
planation to the nature of this fluid is found in the standard model 
of cosmology ACDM, where the role of the DE is played by a cos- 
mological constant A obeying the equation of state pa = — pa- 
Although perfectly consistent with all the aforementioned observa- 
tions, ACDM still lacks of appeal from a purely theoretical point of 
view. In fact, if we belive the cosmological constant to be the zero 
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point energy of some fundamental quantum field, its introduction 
in the Friedmann equations requires a fine-tuning of several tens of 
orders of magnitude (depending on the energy scale we choose to 
be fundamental in our theory) spoiling the naturalness of the whole 
ACDM picture. 

Another issue we encounter when dealing with the standard cos- 
mological model is the so called coincidence problem, that is, the 
difficulty to explain in a natural way the fact that today's matter 
and dark energy densities have a comparable value although they 
evolved in a completely different manner throughout most of the 
history of the universe. 

In an attempt to overco me these two difficulties of ACDM, 
IJimenez & Marital d2009h introduced the Vector Dark Energy 
(VDE) model, where a cosmic vector field plays the role of a dy- 
namical dark energy component, replacing the cosmological con- 
stant A. Besides being compatible with supernovae observations 
and CMB precision measurements, this scenario has the same num- 
ber of free parameters as ACDM. Moreover, the initial value of 
the vector field (which is of the order of 10~ 4 A/J3, a scale that 
could arise naturally in inflation) and its global dynamics ensure 
the model to overcome the standard model's naturarlness problems. 
In the present work we study the impact of this VDE model on 
structure formation and evolution by means of a series of cosmo- 
logical TV-body simulations, analyzing the effects of this alterna- 
tive cosmology in the deeply non-linear regime and highlighting its 
imprints on cosmic structures, in particular, emphasizing the differ- 
ences emerging with respect to the standard model ACDM. To be 
able to disentangle the effects due to the different parameters from 
those induced by the different dynamics of the background in VDE 
cosmology, we have also run a set of simulations for a ACDM-vde 
cosmology, i.e. a model that embeds VDE parameters (presented in 
Section^ in a standard ACDM picture. 

Although ruled out by current cosmological constraints, this 
model provides nonetheless an interesting case study that allows us 
to shed light on the effects of these two cosmological parameters 
on structure formation in the VDE model. 

The paper is organized as follows. In Section|2]we briefly in- 
troduce the VDE model, discussing its most important mathemati- 
cal and physical characteristics. In Section|3]we describe the setup 
as well as the modifications to the code and the initial conditions 
necessary to run the A-body simulation. In the two Sections|4]and 
[5]we will present a detailed analysis of the results, focusing on the 
main differences of the VDE models to the standard ACDM cos- 
mology, first analyzing the large-scale structure and then (cross- 
comparing properties of dark matter haloes. A short summary of 
the results obtained and a discussion on their implications is then 
presented in section Section|6] 



2 THE MODEL 

In this section we will provide the basic mathematical and physical 
description of the VDE model. For more details and an in-depth 
discussion on the results obtained and their derivation we refer to 
iBeltran Jime nez & Maroto (2008). The action of the proposed vec- 



M p being the Planck mass. 



tor dark energy model can be written as: 
R 1 



S 



-F F^ 



1 



(1) 



where /f M „ is the Ricci tensor, R = the scalar curvature 

and F ^ — d^Av — dvA^. This action can be interpreted as a 
Maxwell term for a vector field supplemented with a gauge-fixing 
term and an effective mass provided by the Ricci tensor. It is in- 
teresting to note that the vector sector has no free parameters nor 
potential terms, being G the only dimensional constant of the the- 
ory. This is one of the main differences of this model with respect 
to those based on scalar fields, which need the presence of potential 
terms to be able to lead to late-time accelerated expansion. 

The classical equations of motion derived from the action (0 
are the Einstein and vector field equations given by: 



OA ii + R^A V 



0. 



(2) 
(3) 



where T ' M „ is the conserved energy-momentum tensor for matter 
and radiation (and/or other possible components present in the Uni- 
verse) and Tz„ is the energy-momentum tensor coming from the 
vector field sector (and that is also covariantly conserved). In the 
following, we shall solve the equations of the vector field during 
the radiation and matter eras, in which the contribution of dark en- 
ergy is supposed to be negligible. In those epochs, the geometry 
of the universe is well-described by the flat Friedmann-Lemaitre- 
Robertson-Walker metric: 



ds 2 = dt 2 — a(t) 2 dx 2 



(4) 



For the homogeneous vector field we shall assume, without lack of 
generality, the form A M = (Ao(t),0,0,A z (t)) so that the corre- 
sponding equations read: 



A + 3HA - 3 [2H 2 + H ) A = 



A z + HA Z -2(H + 3H 2 ) A z = 0, 



(5) 
(6) 



where H — a/a. These equations can be easily solved for a power 
law expansion with H = p/t, in which case we obtain the follow- 
ing solutions: 



A (t) = A+t a + +Aot a ~, 
A z {t) = Att 2p +A z t 1 - Sp . 



(7) 
(8) 



with q± = |(1 - 3p ± \/33p 2 - 18p + 1) and A± and Af 
constants of integration. Thus, in the radiation dominated epoch 
(p = 1/2) we have the growing modes Aq ^constant and A z oc t, 
whereas for the matter dominated epoch (p = 2/3) we have 
A oc f(- 3 +^ / 33)/ 6 an d A z oc t 4/3 . Concerning the energy den- 
sities, the corresponding expressions are given by: 



PA = \h 2 A 2 + 3HA A - iig, 



PA Z = ^ (4H 2 A 2 Z - AHA Z A Z + A 2 



(9) 
(10) 



At this point, it is interesting to note that when we insert the full 
solution for A z given in <[S]> in its corresponding energy density, we 
obtain: 



PA Z = 



2a s 



(25p 2 - Wp + 1) 



(11) 
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Figure 1. Left plot: Evolution of the energy densities. Dashed (red) for radiation, dotted (green) for matter and solid (blue) for vector dark energy. We also 
show for comparison the cosmological constant energy density in dashed-dotted line. We see the scaling behaviour of the cosmic vector in the early universe 
and the rapid growth of its energy density contribution at late times when approaching the final singularity. Right plot: Evolution of the temporal component 
of the vector field where we see that it takes a constant value at very high redshifts so that the cosmological evolutions is insensitive to the precise redshift at 
which we set the initial value of the cosmic vector. 



so that the mode A f does not contribute to the energy density. That 
way, even though A z grows with respect to Aq, the correspond- 
ing physical quantity, i.e. its energy density, decays with respect 
to that of the temporal component. It is easy to check that the ra- 
tio Pa z /pa decays as a~ 4 in the radiation era and as a~ 6 ' 37 in 
the matter era so that the energy density of the vector field becomes 
dominated by the contribution of the temporal component. That jus- 
tifies to neglect the spatial components and deal uniquely with the 
temporal one. 

On the other hand, the potential large scale anisotropy gener- 
ated by the presence of spatial components of the vector field is 
determined by the relative difference of pressures in different di- 
rections p|| and p±, that is given by: 



P\\ ~P± 



4HA Z A Z + Ai 



(12) 



This expression happens to be proportional to pA z so that we have 
that (p|i — p±)/pa will decay as the universe expands in the same 
manner as pa z and the large scale isotropy of the universe sug- 
gested by CMB observations is not spoiled. Hence, in the following 
we shall neglect the spatial components of the vector field and we 
shall uniquely consider the temporal one, since it gives the dom- 
inant contribution to the energy-momentum tensor of the vector 
field. However, we should emphasize here that this does not re- 
sult in effectively having a scalar field. As commented before, for a 
minimally-coupled scalar field, one needs to introduce a certain po- 
tential (that will depend on some dimensional parameters) to have 
accelerated expansion, whereas in the VDE model, we get acceler- 
ated solutions with only kinetic terms and without introducing any 
new dimensionful parameter. 

The energy density of the vector field is given by: 



PA = PAo(l + z) 



(13) 



with k — 4 in the radiation era and k = (9 — v33) /2 ~ —1.63 in 
the matter era. We can also calculate the effective equation of state 
for dark energy as: 



W DE 



PA 
PA 



-3 (§H 2 + ftf) Al + HA A 



2 A 2 
2 A) 



H 2 Al + 3HA A - \A\ 



(14) 



Again, using the approximate solutions in ((7), we obtain; 



W D E = 



3V33-13 
\/33-15 



radiation era 
-0.457 matter era 



(15) 



From the evolution of the energy density of the vector field we see 
that it scales as radiation at early times, so that pa/ pa = const. 
However, when the Universe enters its matter era, pa starts grow- 
ing relative to pm eventually overcoming it at some point, in which 
the dark energy vector field would become the dominant compo- 
nent. From that point on, we cannot obtain analytic solutions to the 
field equations and we need to numerically solve the corresponding 
equations. In Fig.[T]we show such a numerical solution to the exact 
equations, which confirms our analytical estimates in the radiation 
and matter eras. Notice that, since Ao is constant during the radia- 
tion era, the solutions do not depend on the precise time at which 
we specify the initial conditions as long as we set them well inside 
the radiation epoch. Thus, once the present value of the Hubble pa- 
rameter Ho and the constant Aq during radiation (which indirectly 
fixes the total matter density £Im) are specified, the model is com- 
pletely determined. In other words, this model contains the same 
number of parameters as ACDM, i.e. the minimum number of pa- 
rameters of any cosmological model with dark energy. 

The vector dark energy model does not only have the min- 
imum required number of parameters, but it also allows to alle- 
viate the so-called naturalness or coincidence problems that most 
dark energy models have. This is so because the required value 
for the constant value that the vector field takes in the early uni- 
verse happens to be ~ 10 _4 Af p . This value, in addition to be rela- 
tively close to the Planck scale, could naturally arise from quantum 
fluctuations during inflation, for instance. On the other hand, the 
fact that the energy density of the vector field scales as radiation 
in the early universe also goes in the right direction of alleviat- 
ing the aforementioned problems because the fraction of dark en- 
ergy during that period remains constant. Moreover, said fraction 
is Q, A arly = pa/ Pr — 1CP 6 , which, again, is in agreement with 
the usual magnitude of the quantum fluctuations produced during 
inflation. 

After dark energy starts dominating, the equation of state 
abruptly falls towards wde — > — oo as the Universe approaches 
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a finite time t en( j. As shown in Fig.[2] during the cosmological evo- 
lution the equation of state crosses the so-called phantom divide 
line, so that we have wde{z — 0) < — 1. The final stage of the 
universe in this model is a singularity usually called Type III or Big 
Freeze, in which the scale factor remains finite, but the Hubble ex- 
pansion rate, the energy density and the pressure diverge. This is 
a distinct feature of the VDE model as compared to quintessence 
fields for which the equation of state is restricted to be > — 1 so that 
no crossing of the phantom divide line is possible. In fact, for a dark 
energy model based on scalar fields, one needs either non-standard 
kinetic terms involving higher derivative terms in the action or the 
presence of several interacting scalar fields to achieve a transition 
from w > —1 to a phantom behaviour (w < —1). In either case, 
non-linear derivative interactions or multiple scalar field scenarios, 
additional degrees of freedom are introduced, whereas the VDE 
model is able to obtained the mentioned transition with only the 
degree of freedom given by the temporal component of the vector 
field. 

Notice that in the VDE model the present value of the equa- 
tion of state parameter Wo = —3.53 is radically different from that 
of a cosmological constant (cf. Fig. 1, where the redshift evolution 
of w(z) is shown in the range of our simulations). The values of 
other cosmological parameters also differ importantly from those 
of ACDM (see Table[j]>. Despite this fact, VDE is able to simulta- 
neously fit supernovae and CMB dat a with comparable goodness to 
ACDM ( iBeltran Jimenez & Marotol J2008t) . iBeitran Jimenez et all 
J2009h ). 

In particular, for CMB0, the \ 2 f° r the best fit parameters for 
ACDM is 48.3, wheres for the VDE model we obtain \ 2 = 51-8 
for the parameters used to run our simulations. Thus, even though 
the equation of state evolution is as the one shown in Fig. [2] the 
VDE model provides good fits to observations, as shown in Fig. [3] 

This might seem to be surprising if we notice that the present 
equation of state for the VDE model is u>o — —3.53, which is far 
from the usual constraints on this parameter obtained from cosmo- 
logical observations. Such constraints are usually obtained by as- 
suming a certain parametrization for the time variation of the dark 
energy equation of state. However, the different used parameteri- 
zations are normally such that ACDM is included in the parameter 
space. If we look at Fig. [2] we can see that the evolution of the 
equation of state for the VDE model crucially differs from those 
of ACDM or quintessence models and, indeed, it cannot be prop- 
erly described by the most popular parameterizations. This means 
that we cannot directly apply the existing constraints to the VDE 
model, but a direct comparison of its predictions to observations is 
required. 

As a final remark, in the simulations we will not include inho- 
mogeneous perturbations of the vector field, but only the effects of 
having a different background expansion will be considered. 

In Fig. [4] we show the matter power spectrum for both ACDM 
and the VDE model. The differences can be ascribed to the fact of 
having different cosmological parameters that change the normal- 
ization and the matter-radiation equality scale k eq , which are the 
only two differences observed. Notice that the transfer function is 
the same in both cases, since the slopes before and after the k eq 
are the same, so that we do not expect strong effects at early times 
which could affect the evolution of density parameters. 



We use the binned data of WMAP7 
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Figure 2. Evolution of dark energy equation of state where we can see the 
crossing of the phantom divide line and its evolution towards — oo as the 
final singularity is approached. 




Figure 3. ACDM and VDE best fit values for the CMB spectrum versus 
WMAP7 data, showing the viability of VDE cosmology. The same curve is 
also shown for a ACDM-vde cosmology, which is clearly ruled out. 



3 THE TV-BODY SIMULATIONS 

In this section we will explain the (numerical) methods used in this 
work, with a particular emphasis on the necessary modifications of 
the standard TV-body and halo finding algorithms, also describing 
the procedures followed to test their accuracy and reliability. 



3.1 Set-up 

The TV-body simulations presented in this work have been car- 
ried out usi ng a suitably m odified version of the Tree-PM code 
GADGET-2 l lSpringelll2005h . It has been also necessary to generate 
a particular set of initial conditions to consistently account for the 
VDE-induced modifications to the standard paradigm. In Table Q] 
we show the most relevant cosmological parameters used in the 
different simulations. For the VDE model, we have used the value 
of Q.m provided by the best fit to SNIa data, whereas the remaining 
cosmological parameters have been obtained by a fit of the model 
to the WMAP7 da taset. For ACDM we used the Multidark Simu- 
lation dPrada etaL 



201 1) cosmological parameters with a WMAP7 



as normalization ( Lar son et al 1l20ll 



In addition, we also simulated a so-called ACDM-vde model, 
which implements the VDE values for the total matter density Qm 
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Table 1. TV-body settings and cosmological parameters used for the 
GADGET-2 simulations, the two 500h, -1 Mpc and the two lh _1 Gpc have 
the same initial random seed (in order to allow for a direct comparison of 
the halo properties) and starting redshift z s tart = 60. The number of par- 
ticles in each ran was fixed at 512 3 . The box size B is given in units of 
h~ 1 Mpc and the particle mass in 10 11 h~ 1 M.Q. 



Simulation 






<T8 


h 


B 


in,, 


2XVDE-0.5 


0.388 


0.612 


0.83 


0.62 


500 


1.00 


2XVDE-1 


0.388 


0.612 


0.83 


0.62 


1000 


8.02 


ACDM-0.5 


0.27 


0.73 


0.8 


0.7 


500 


0.69 


ACDM-1 


0.27 


0.73 


0.8 


0.7 


1000 


5.55 


ACDM-0.5vde 


0.388 


0.612 


0.83 


0.7 


500 


1.00 


ACDM-1 vde 


0.388 


0.612 


0.83 


0.7 


1000 


8.02 



and fluctuation amplitude as in an otherwise standard ACDM pic- 
ture. 

In particular, we want to be able to determine the impact of 
the different parameters on cosmological scales, with a particular 
emphasis on the very large structures and the most massive clus- 
ters, where observations are sta rting to clash with the predictions of 
the cu rr ent standard model (see Jee et al. l d2009h . lBaidi & Pettorinol 



2011 



), lHovleetai] J201 lh . lCarlesi et alJoOllI) and lEnqvist et al 



20111) ). Therefore, we need to determine whether the results de- 
rived from our VDE simulations can be solely attributed to its ex- 
tremely different values for the cosmological parameters or actu- 
ally by the presence of the cosmic vector field. In other words, we 
want to separate the signatures of the dynamics-driven effects from 
the parameter-driven ones, with a focus on large scale structures, 
where the imprints are stronger and more clearly connected to the 
cosmological model. We chose to run a total of eight 512 3 particle 
simulations summarized in Table[T]and explained below: 

• two VDE simulations, i.e. a 500 ft -1 Mpc and a 1 h~ 1 Gpc 
box, 

• two ACDM simulations with the same box sizes and initial 
seeds as the VDE runs above, 

• two more VDE simulations with a different random seed, 
again one in a 500 h~ Mpc and one in a 1 /i _1 Gpc box (both 
serving as a check for the influence of cosmic variance), and 

• two ACDM-vde simulations, one again in a 500 h' 1 Mpc and 
one in a 1 h~ 1 Gpc box. 

All runs were performed on 64 CPUs using the MareNostrum 
cluster at the Barcelona Supercomputing Center. Most of the results 
we will discuss and analyze here are based on the 500 h" 1 Mpc 
simulations as they have the better mass resolution. The 1 h~ 1 Gpc 
runs primarily serve as a confirmation of t he results and have al- 
ready been discussed in lCarlesi et al.l d201 ll) . respectively. 
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Figure 4. Linear matter overdensity power spectra at z = and z = 60 for 
VDE, ACDM and ACDM-vde plotted versus wavenumber k. Vertical solid 
thick black lines refer to the fc-space interval covered by the 500 h Mpc 
simulations whereas the thin ones refer to the 1 h~ 1 Gpc one. All matter 
power spectra at z = have been normalized to the erg values shown in 
Table[T]and then rescaled to z = 60 via the linear growth factor. We notice 
that for /« < 0.05hMpc _1 , ACDM andACDM-vde have more power than 
VDE, whereas on smaller scales the opposite is true. We also note that due 
to the different value of erg normalization the ACDM-vde P(k) is slightly 
larger than the ACDM one at z = while the different growth factor, 
which is larger in the ACDM-vde cosmology, affect the setting of the initial 
conditions, where the latter power spectrum lies below the former. 
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Figure 5. The ratio of the Hubble function H^h,- 1 for VDE and ACDM- 
vde to the standard ACDM one. At earlier times VDE undergoes a rela- 
tively faster expansion compared to ACDM, whereas the opposite is true at 
smaller zs. On the other hand, ACDM-vde cosmology is characterized by a 
slower relative expansion throughout the whole history of the universe. 



3.2 Code Modifications 

In the following paragraph we are going to describe the procedures 
followed to implement the modifications needed in order to run our 
TV-body simulations consistently and reliably. This is in principle a 
non-trivial issue, since, as described in Sectionf2] we need to incor- 
porate a large number of different features that affect both the code 
used for the simulations and the initial conditions. 

In particular, we have to handle with care three features that 
distinguish it from ACDM, i.e.: 

• the matter power spectrum P(k,z) (shown in Fig. [4]l and its 



normalization as, 

• the expansion history H(z) (see Fig.[5J, and 

• the linear growth factor D + (z) (cf. Fig.[6]l. 

Whereas the first and last point affect the system's initial con- 
ditions, the second one enters directly into the TV-body time inte- 
gration, and has to be taken into account by a modification of the 
simulation code. 
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3.2.1 Initial conditions 

To consistently generate the initial conditions for our simulation, 
first we normalized the perturbation power spectrum depicted in 
Fig.|4]to the chosen value for a$ at z = 0. Therefore, we normal- 
ized VDE and ACDM-vde inital conditions to as = 0.83 while for 
ACDM we used the WMAP7 value as = 0.8. Using the respective 
linear growth factors, we rescaled the P(k) to the initial redshift 
z = 60 where then the particles' i nitial velocities a nd positions 
were computed using the Zeldovich jZeldovich|[l970h approxima- 
tion. 

We emphasize here that the main goal of our analysis is to find 
and highlight the main differences of the VDE picture with respect 
to the standard one: therefore, the choice of these different nor- 
malization parameters has to be understood as unavoidable as long 
as we want the models under investigation to be WMAP7 viable 
ones. Needless to say, in this regard the ACDM-vde cosmology 
must be considered only as a tool to disentangle parameter-driven 
effects from the dynamical ones, not being a concurrent cosmolog- 
ical paradigm we want to compare VDE to. 
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Figure 6. Ratio of the growth function to the expansion factor D(a)/a as 
obtained from the 500 h^ 1 Mpc box simulations versus the analytical one. 
The results shown an agreement between the theoretical expectation and 
the numerically computed one within the 2% level. The results from the 1 
h~ 1 Gpc box simulations are not shown since they perfectly overlap with 
the ones presented here. 



3.4 Halo Finding 



3.2.2 Hubble expansion 

As pointed out bv lLi & Barrowl ( l201 ll) . the expansion history of the 
universe has a very deep impact on structure formation and in par- 
ticular the results of an A/-body simulation, as it affects directly 
every single particle through the equations of motion written in 
comoving coordinates. In Fig. [5] the ratios of the Hubble expan- 
sion factors for VDE and ACDM-vde to the standard ACDM value 
are shown; we see that different models are characterized by dif- 
ferences up to the 20% in the expansion rate. To implement this 
modification, we replaced the standard computation of H(a) in 
GADGET- 2 with a routine that reads and interpolates from a pre- 
computed table. 



3.3 Code Testing 

To check the reliability of the modifications introduced into the 
simulation code and during the generation of the initial conditions, 
we have confronted the theoretica l linear growth fac tor, computed 
using the Boltzmann-code CAMB l lLewis et alfeOOOh with the ones 
derived directly from the simulations. 

As shown in Fig. [6] our results yield an agreement within the 
1% level, which proofs the correctness of our modifications as well 
as illustrating (again) the differences in structure growth between 
the models. 

We would like to note that for consistency reasons, both when 
calculating the CAMB and the numerical value for the growth factor, 
we have used the expression 



D + (z) 



P{z,k ) 
P{zo, ko) 



(16) 



In order to identify haloes in our simulation we have run the open 
so urce MPI+OpenMP hybrid halo finder AHE0 described in detail 
in lKnollmann & Knebel 120091 ). AHF is an improvement of the MHF 
halo finder dGill et al.ll2004l) and has been extensive ly compared 
agains t practically all other halo finding methods in iKnebe et al .1 
d201 lh . AHF locates local overdensities in an adaptively smoothed 
density field as prospective halo centres. For each of these den- 
sity peaks the gravitationally bound particles are determined. Only 
peaks with at least 20 bound particles are considered as haloes and 
retained for further analysis. 

But the determination of the mass requires a bit more elabora- 
tion as it is computed via the equation 



M(R) = A x p c (z) x ~^-R 



(17) 



where we applied A = 200 as the overdensity threshold and p c (z) 
refers to the critical density of the universe at redshift z. In this 
way M(R) is defined as the total mass contained within a radius 
R, corresponding to the point where the halo matter density p(r) is 
A times the critical value p c . Using this relation, particular care has 
to be taken when considering the definition of the critical density 



Pc(z) = 



3H 2 (z) 



(18) 



8ttG 

because it involves the Hubble parameter, that differs substantially 
at all redshifts in the two models. This means that, identifying the 
halo masses, we have to take into account the fact that the value of 
pc(z) changes from ACDM to VDE. This has been incorporated 
into and taken care of in the latest version of AHF where Hvde{z) 
is being read in from a precomputed table, too. 

We finally need to mention that we checked that the objects 
obtained by this (virial) definition can be compared across differ- 
ent cosmological models and using different mass definitions. To 



where ko is a fixed scale whithin the linear regime and zo is the 
initial redshift of the simulation. 



3 AHF stands for Amiga Halo Finder, to be downloaded freely from 

http : / /www . popia . f t . uam . es/AMIGA 
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this extent we studied the ratio between two times kinetic over po- 
tential energy r\ = 2T/\U\ confirming that at each redshift under 
investigation here the distributions of r\ in ACDM and VDE are 
actually comparable (not presented here though), meaning that the 
degree of virialisation (which should be guaranteed by Eq. (1171 )) is 
in fact similar. We therefore conclude that our adopted method to 
define halo mass (and edge) in the VDE model leads to unbiased 
results and yields objects in the same state of equilibrium as is the 
case for the ACDM haloes. Please note that this test does not guar- 
antee that all our objects are in fact virialized; it merely assures us 
that the degree of virialisation is equivalent. We will come back to 
this issue later when selecting only equilibrated objects. 



4 LARGE SCALE STRUCTURE AND GLOBAL 
PROPERTIES 

In the following section, we will discuss the global properties of 
large scale structures identified in our simulations. Using all of our 
sets of simulations for ACDM, ACDM-vde and VDE we will dis- 
entangle parameter-driven effects from those due to the different 
dynamics of the background expansion, which uniquely charac- 
terize VDE and therefore are worth pointing out in the process of 
model selection. 



4.1 Density Distribution 

In Fig. [7] we show the colour coded density field for the particle 
distribution at redshift z = 0, for a 120 x 120 h' 1 Mpc 2 slice at 
the box center for the three 500 h^ 1 Mpc simulations projected on 
the x-z plane. As expected, we observe that the most massive struc- 
tures' spatial positions match in the three simulations, although in 
the ACDM-vde and VDE a large overabundance of objects with 
respect to ACDM, as we could expect due to the higher fljj. This 
observation will be confirmed on more quantitative grounds in the 
analysis earned in the following sections, especially when refering 
to the study of the cumulative mass function. 

4.2 Matter Power Spectrum 

In Fig.[8]we show the dark matter power spectrum P(k) at redshifts 
z = 0, 1, 3, 4 computed for the VDE-0.5, ACDM-0.5 and ACDM- 
vde-0.5 simulations. For clarity, we do not show the 1 h~ Gpc sim- 
ulations; however, we have checked their consistency with the 500 
h^ 1 Mpc runs. We note that at all redshifts the differences already 
seen in the input power spectra are preserved (cf. Fig. @J, mean- 
ing that the VDE model has less power than ACDM on the large 
scales, whereas the opposite is true for small scale. This particu- 
lar shape of the P(k) is a peculiar fea ture of VDE cosm ology, as 
other kin ds of dynamical qu intessence jAlimi et"al J2010h and cou- 
pled DE tealdi et al]|2010h show completely different properties; 
with less power (in the former case) or a ACDM-type of behaviour 
(in the latter) on small scales. At higher and intermediate redshifts 
ACDM-vde shows almost no differences from ACDM, as expected 
since the former is normalized to a lower initial value with respect 
to the latter and therefore needs to equal it before eventually over- 
coming it at smaller z's, as imposed by the larger erg normalization. 
The effects of the different growth factor in this model start to be- 
come evident only at z < 1, where we see that the ratio of the 
P(k) starts to increase. Whereas the ratio of VDE to ACDM for 
k < 0.05 h h~ x Mpc -1 is substantially unaltered at all redshifts, 



small scales are affected by non-linear effects, eventually distorting 
its shape. 

4.3 Halo Abundance 

In the following subsection we will study the abundance of massive 
objects at different redshifts. Highlighting the differences arising 
among the three models in the different mass ranges, we want to 
study VDE's peculiar predictions for the massive cluster distribu- 
tion and highlight its distinction from ACDM. 

To this extent, we compare in figure Fig.QT]the three different 
mass functions at z — computed for the lfe -1 Gpc simulations 
versu s the ones obtained using the Tinker formula dTinker et al.l 
l2008h . In the ACDM and ACDM-vde cases the theoretical val- 
ues are obtained using the standard parametrization, whereas in the 
VDE case we performe a non linear fit of the mass function to the 
numerical data. We find that in the VDE case the four parameters 
of the Tinker mass functions take the values (assuming the stan- 
dard notation) A = 0.105, a = 0.75, b = 4.42 and c = 1.48; thus 
differing substantially from the usual ones (which are A = 0.187, 
a — 1.47, b = 2.57 and c = 1.19). This difference is most likely 
due to a different value for the linear critical overdensity parameter 
S c ; however, a complete understanding of this aspect would require 
a deeper knowledge of the linear perturbation theory in VDE cos- 
mology, which is still under investigation. 

In Fig.[9]we show the (cumulative) mass functions for the three 
models at z = 0, 1, 2, 4, as computed from the VDE-0.5, ACDM- 
0.5 and ACDM-0.5-vde simulations; the correspon ding VDE-1, 
ACDM -1 and ACDM-vde- 1 results can be found in ICarlesi et al.l 
d201 ll) ; they are not shown here again as they do not provide any 
new insights and rather confirm (and extend) the results to be drawn 
from the 500/i _1 Mpc boxes, respectively: We note that the VDE 
cosmology is characterized by a larger number of objects at all the 
mass scales and redshifts, outnumbering ACDM by a factor con- 
stantly larger than 2. In particular, this enhancement can be seen 
for the very large masses, where at low z the VDE/ ACDM ratio 
reaches values of ~ 10. Although this value of the ratio seems 
to be a mere result of the cosmic variance, due to the low num- 
ber of haloes found in this mass range, the computation of the 
mass function for the second 500 h~ x Mpc VDE realization and 
the l/i _1 Gpc simulations makes us believe that the expected en- 
hancement in this region must be at least a factor 5. 

Interestingly enough, ACDM-vde has comparable character- 
istics to VDE, which leads to the conclusion that the substantial en- 
hancement in structure formation is mainly parameter-driven, i.e. 
due to the overabundance of matter and higher normalisation of 
matter density perturbations. Although this first observation may 
seem in contrast with what we have found in Section l4~2l where 
we have noticed that VDE has less power on large scales in com- 
parison to ACDM, we have to take into account that, in the hi- 
erarchical picture of structure formation, objects on small scales 
form first to subsequently give birth to larger ones. This means, in 
our case, that more power for large fc-values should be regarded 
as an important source of the overall enhancement together with 
the overabundance of matter, as already pointed out in the previous 
discussion. The evolution of the mass functions at different redshift 
allows us to disentangle the effect of the modified expansion rate; 
at higher redshift, in fact, both the ACDM and ACDM-vde mass 
functions are suppressed with respect to the VDE model, mostly 
because of the lack of power on small scales. These stronger initial 
fluctuations eventually trigger the earlier start of structure forma- 
tion, but - as time passes - the effect of the increased expansion 
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Figure 7. Projected density for ACDM,ACDM-vde and VDE showing a 120 X 120 h~ 2 Mpc 2 slice at the box center in the 500 /i -1 Mpc box at z = 
projected on the x-z plane. Bright areas are associated with matter whereas underdense regions are denoted by darker, black spots in the projected box. Results 
for the VDE-1, ACDM-1 and ACDM-vde-1 simulations are not shown since the colour coding does not provide useful insights on the different clustering 
patterns on smaller scales. 



rate shown in Fig. [5]for the VDE cosmology suppresses structure 
growth, leading to a mass function below the ACDM-vde curve at 
around redshift one. At this point, the VDE expansion rate starts 
decreasing with respect to the ACDM one, comparatively enhanc- 
ing very large structure growth and eventually causing the two mass 
functions to be (nearly) indistinguishable at z = 0. In an attempt 
to disentangle further the effects of the new Qm to those due to the 
higher normalization of the matter power spectrum, we searched for 
a redshift zq at which ACDM and VDE have an identical erg value, 
which would allow us to single out the impact of the different CIm 
parameter alone. Integrating numerically the power spectra of the 
different simulations' snapshots, we have found zq — 0.18, where 
erg is 0.65 1 for both VDE and ACDM. Analyzing the power spectra 
and mass functions at this redshift we can conclude that the main 
effects of VDE, such as the overabundance of objects, are only due 
to the larger Qm, being the impact of the different normalization of 
the matter fluctuations practically negligible. 

Furthermore, if we look at Fig. [10] where we show the evo- 
lution with redshift of the number density of objects above the 
M = 10 14 /i _1 M o threshold, we observe that the most massive 
structures in the two cosmologies form at comparable rates. This 
seems to suggest that in the VDE picture there is a subtle balance 
between the formation of new small haloes and their merging into 
more massive structures. Such an effect comes as no surprise if 
we again take into account that this model has two main opposite, 
different features that affect the formation of structures: a strong 
suppression on all scales induced by the faster expansion of the 
universe for a large redshift interval and an enhancement due to a 
higher density of matter and a larger power on the small scales. 

An interesting consequence of this kind of behaviour is that 
the VDE overabundance of massive objects may address some re- 
cent observational tensions of ACDM; namely, the high redshift 
of reionization and the presence of extremely massive clusters at 
2 > 1. Recent microwave background observations seem to prefer 
a high reionization redshift, around z m 10 combined with a lower 
normalization of the matter perturbation s, gg ~ 0.8; whereas sim- 
ulations have shown (see, for example, iRaicevic et alj|20llh that 
early reionization can be achieved only for erg = 0.9 or larger. In 
VDE, the appearence of dark matter haloes with masses larger than 



1O 12 /i -1 M0 as early as z = 7 (while equivalent structures appear 
in ACDM only for z > 5) might imply also a larger z re ion, pro- 
vided the hierarchycal picture of structure formation holds also in 
VDE at smaller mass scales. On the other hand, the exi stence of 
M > 5 X 10 14 fe~ 1 Mp ) clu s ters at z > 1 (as d iscussed in |jee et all 



120091 ; iBrodwin et a"i]|20irj| ; iFolev et alj|201lh has also been con- 
sidered by many authors (e.g., Hovle et al .l201lHBaldi&Pettorinol 



1201 ll ; lBaldil2012l ; lEnavist et alj|201 ll) as a serious challenge to the 
standard ACDM paradigm; for a more thorough discussion of this 
issue in the conte xt of VDE cosmolog y we refer to aforementioned 
articles as well as ICarlesi et al]j2oT lj). However, the comparison to 
the ACDM-vde paradigm also shown in Fig.llOlshows that indeed 
VDE acts as a source of suppression of structure growth with re- 
spect to the enhancement triggered by the increase in erg and Qm. 
This effect is i ndeed a general result of uncoupled dyn amical dark 
energy models dGrossi & Springell2009tlLi et al.l201 lj) as the pres- 
ence of a larger fraction of dark energy at high 2 enhances Hubble 
expansion (as shown in Fig. [5} preventing a stronger clustering to 
take place. 

In our case, it is also important to point out that the overpre- 
diction of objects at z = may represent a shortcoming of the 
model, as observations on the cluster num ber mass function seem 
to be in contrast w ith s uch a prediction (seelVikhlinin et all d2009T ). 
IWen et al.l d2010h and lBurenin & Vikhlinir] d2012h l. Furthermore, 
we have to keep in mind that these results assume a ACDM fidu- 
cial model, while the use of a different cosmology requires a careful 
handling of the data and does not allow a straightforward compar- 
ison to the observations, as they are affected by model-dependent 
quantities like comoving volumes and mass-temperature relations. 



4.4 Void Function 

In order to identify voids, our voidfinder starts with a selection of 
point like objects in 3D. These objects can be haloes above a certain 
mass or a certain circular velocity or galaxies above a certain lumi- 
nosity. Thus the detected voids are characterized by this threshold 
mass, circular v elocity or luminosity . Other voidfinders use differ- 
ent approaches dColberg et al.ll2008h . The void finding algorithm 
does not take into account periodic boundary conditions used in 
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Figure 8. Power Spectra at redshifts z = 0, 1, 2, 4; solid lines are for for ACDM-0.5, dotted for ACDM-0.5-vde and dashed for VDE-0.5. The results from the 
1 /i _1 Gpc simulations are not shown as they simply overlap to the present ones on the smaller-fc end, without providing further insights on the small scales, 
where we expect non linear effects to dominate. 




Figure 10. Number density evolution for objects more massive than 10 
h~ 1 VIq?s a function of redshift. The larger amount of massive clusters at 
higher redshift is a distinctive feature of VDE cosmology. 



Figure 1 1. The numerical halo mass function for the three 1 h ~ 1 Gpc simu- 
lations at z = 0, showing the numerical results versus the Tinker theoretical 
values. 
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Figure 9. Mass function for ACDM(solid lines), VDE (dashed lines) and ACDM-vde (dotted lines) models at different redshifts, computed for the 500 
h Mpc box simulations. We have also verified that the corresponding values computed for the 1 /i _1 Gpc simulations overlap to the ones shown here for 
M > 10 13 /i — 'Mq, except for a smoother high-mass end. In the lower panels of the plots, VDE and ACDM-vde to ACDM ratios are represented by dotted 
lines while VDE to ACDM-vde are shown using dash-dotted lines. 



numerical simulations. Therefore, we have periodically extended 
the simulation box by 50 h^ 1 Mpc. In this extended box we repre- 
sent all haloes with a mass above the threshold of 5 x 10 12 /i^ 1 Mq 
as a point. In this point distribution we search at first the largest 
empty sphere which is completely inside the box. To find the other 
voids we repeat this procedure however taking into account the pre- 
viously found voids. We allow that newly detected voids intersect 
with previously detected ones up to 25% of the radius of a previ- 
ously detected larger void. 



In Fig.[T2]we show the cumulative number of voids with radius 
larger than i? vo id the center of which is in the original box. One 
can clearly see that for a given void radius there exist more voids 
in the ACDM than in th ACDM-vde and VDE models. The void 
distribution reflects the behaviour of the mass function shown in 
Fig. [9] At redshift z — there exist less haloes with mj, > 5 x 
1O 12 /i _1 M in the ACDM model than in the other two models. 
Thus on average larger voids are expected. 



4.5 Growth Index 

The growth of the perturbations can be related to the evolution of 
the matter density parameter by the general relation 



7 ( Q ) _ dln(<5(q)) 
dln(a) 



(19) 



In the standard ACDM cosmology, the exponent 7(a) can be 
approximated by a constant value 7 ~ 0.55, although a more de- 
tailed calculat ion shows that this numbe r is actually redshift de- 
pendend (see iBueno Belloso et al]|201 lh . In Fig.|T3]we show the 
evolution of this growth index 7(2) computed from our VDE-0.5, 
ACDM-vde and ACDM-0.5 simulations. As expected, we do ob- 
serve that in VDE structure formation is generally suppressed with 
respect to ACDM as an effect of the faster expansion rate. This 
statement is true until z ~ 1.5, when the ratio Hvde / Hacdm 
starts decreasing causing the steep increase in the growth index, 
eventually reducing again as soon as vector dark energy enters 
into the phantom regime (see Section [2), undergoing an acceler- 
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Figure 12. Void function for VDE-0.5, ACDM-0.5 and ACDM-vde-0.5 at 
z=0. For the 500/i -1 Mpc box we show the cumulative number of empty 
spheres of radius R which do not contain any object with mass larger than 
5 x lO 12 h- 1 M . 



pendence on z unaltered. Therefore 7(2) can be effectively used 
as a tool for model selection, embodying effectively VDE's pe- 
culiar equation of state us(z) and expansion history. Current ob- 
servational b ounds on 7 constrain only weakly its value at high 
z's (see e.g. iNesseris & Perivolaropou los 2008) or even favour a 
higher 7(2 = 0) 1 Basilako'sf 2012h in contrast to theoretical cal- 
culations based on ACDM. In any case, it will surely be some- 
thing t o be looked at in the near future, when deep surveys like 
Euclid l lLaureiis et al. 20 111) will provid e stringent constraints on 



buclid (Laurens et al.N2Ul 1 | ) will provid e 
this quantity dBueno Belloso etalj|201ll). 



5 DARK MATTER HALOES 

In this section we will discuss properties of (individual) haloes in 
VDE and ACDM. In particular, we will compare the distributions 
of masses, shape parameter, spin parameter, concentrations and for- 
mation redshifts as well as the shape of dark matter density pro- 
files. In this way, we will determine the most important features 
that characterize on the average a single cosmological model. But 
in addition we are also cross correlating haloes in the two models 
studying differences on a one-to-one basis. By this we will be able 
to determine how the properties of a single given structure change 
when switching from one cosmological picture to the other. 
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Figure 13. Growth index in the VDE and ACDM cosmologies from z = 5 
to 2 = 0. Whereas ACDM's growth index has an almost constant behaviour 
with a mild dependence on the redshift, VDE changes dramatically from a 
regime where growth is relatively suppressed (until z sa 1.5) to a relative 
enhancement at earlier times, where 7 becomes larger. 



ated expansion that strongly suppresses structure formation. This 
latter change, which takes place at z ~ 0.5, is reflected by the peak 
of 7(2), which is reached for the same z. Actually, a s stressed by 
different parametrizations dBueno Belloso et al1l201ll) , the growth 
index is extremely sensitive to the value of the equation of state 
u(z), although an explicit form in terms of VDE cosmology still 
has to be found. Indeed, the extremely different behaviour of this 
parameter at different redshifts is an interesting feature that clearly 
distinguishes the two models in a unique way: In fact, parameter- 
induced modification accounts for a ~ 5% change for the value of 
the growth factor, as the comparison among ACDM and ACDM- 
vde suggests. In this case we observe a slight increase of the 
value of 7(2) at all redshifts, due to the increased growth rate in 
ACDM-vde also shown in Fig. [6] However, these changes have 
no impact on the shape of this function, which keeps its mild de- 



5.1 General properties 

To have a reliable description of the general halo properties we 
need to properly select our sample from the catalogues, in or- 
der to include only those objects composed of a number of par- 
ticle sufficient to resolve its int ernal structure without excee ding 
statistical uncertai nty. Following iMunoz-Cuartas et al. I l20ll and 
Prad a et al. we set this number t o approximately 500, even 

though other au thors (see for example iMaccio et al.l d2007l) and 
lBettetal.lj2007|) ) suggest that lower values can be used, too. How- 
ever, since we are dealing with different simulations run with par- 
ticles of different mass, the application of this criterion is not 
straightforward. In fact, since our aim is to compare equivalent 
structures (i.e. structures with the same M200) and not structures 
composed by an identical number of particles we need to choose 
our sample imposing a mass threshold Mth- For the simulations in 
the 500/i -1 Mpc box, we have chosen M th = 5 x W 13 h~ 1 M Q , 
which corresponds to haloes formed by at least 500 particles in 
VDE and ACDM-vde and 715 particles in ACDM; while for the 
larger 1000ft -1 Mpc runs we imposed a M th = 3 x 10 14 /i _1 M o 
limit, i.e. 380 VDE and ACDM-vde particles and 545 ACDM ones. 
In the latter set of simulations, we see that we are including also 
haloes with a ~ 20% less than 500 particles in the VDE and 
ACDM-vde cases; this has been done since in the trade off between 
resolution and sample size we have felt more comfortable using a 
larger number of haloes at the expense of a slight reduction in ac- 
curacy, which will be nonetheless taken into account when analyz- 
ing the results. The total number of haloes that comply with these 
conditions in every simulations, as well as the number of haloes 
that satisfy the relaxation criterion which will be discussed in Sec- 
tion l5.1.3l are shown in Table [2] The state of virialisation of haloes 
will only be taken into account below when investigating the den- 
sity profiles; for the study of the (distributions of the) two-point cor- 
relation functions, the spin, and even the shape of haloes we prefer 
to include even un-relaxed objects as they should clearly stick out 
in the distributions (if present in large quantities). 
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Table 2. Number of haloes above the mass (number) threshold M t h (Nth) 
per simulation. It is also shown the number of relaxed haloes, defined as 
those complying with the criterion introduced in Section |5. 1.31 



Simulation 
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Figure 15. Spin parameter versus the analytical log norm distribution cal- 
culated with the best fit parameters. The fit has been performed using the 
combined sample of haloes above 5 x 10 13 h _1 MQ belonging to the three 
500 /i -1 Mpc boxes with those above the 3 x 10 14 /i _1 M o limit in the 
lft _1 Gpc boxes. 



5.1.2 Spin parameter, shape and triaxiality 

Rotational properties of the haloes can be studied using the so 
called spin parameter A, a dimensionless number tha t measures the 
degree of rotational support of the halo. Following iBullock et al.l 
u200lh . we define it as 



Figure 14. The two point correlation function for objects in the 
500/i — 1 Mpc simulations more massive than 5 X 1O 13 /i _1 M0. 



5.1.1 Correlation function 

To study the clustering properties of the haloes in VDE cosmology, 
we computed the two point correlation function using the defini- 
tion: 



V ^ m{r; Ar) 

JV2 2^ v ( r - Ar) 

i=i v ' 1 



(20) 



where N is the total number of objects above the given mass thresh- 
old in the simulation volume V, and m is the total number of ob- 
jects within a shell of volume v and thickness Ar (of constant log- 
arithmic spacing in r) centered at the ith object. In this case, we 
have limited our analysis to the 500/i -1 Mpc boxes, ignoring the 
l/i _1 Gpc due to their lack of small scale resolution. The results 
are plotted in Fig. [14] where we can see that the £ (r) is slightly 
smaller at all scales in VDE. Although in principle we would ex- 
pect VDE cosmology to have an enhanced clustering pattern due to 
the increased distribution of massive objects observed in the mass 
function, the N~ 2 dependence of the two point correlation func- 
tion drags the total value down, making the final distribution func- 
tion smaller than in ACDM. In fact, a similar behaviour can be ob- 
served for ACDM-vde; with a two point correlation function below 
ACDM at practically all scales. In Table [3] we show the results of 
fitting £(r) to a power law (ro/V) 7 from which we see that VDE is 
characterized by a smaller correlation length ro and a steeper slope 
7- 



A = 



£200 



a/2 M200 V200 -R200 



(21) 



where the quantities L (the total angular momentum), M (to- 
tal mass), V (circular velocity) and R (radius) are all taken 
at the point where the average halo density becomes 200 
times the critical density. Different authors have found (e.g 



times the critical density. Uitt erent autho rs have found (e.g . 
Barnes & Efstathioull987l;IWarren et al.l 19921; ICole & Lacevll99e 



Gardner! l200ll; IBullock et alj l200ll ; iMaccio et al.l l200l boo" : 
Munoz-Cuartas et alj|201ll) that the distribution of this parameter 
is of lognormal type 



P(A) 



Xo-l 



2.7T 



cxp 



ln 2 (A/A ) 



2a 2 



(22) 



even though there are recent claim s that this distribution has to be 
slightly modified jBett et alj|2007l) . 

Fitting the above function to our numerical sample by a non- 
linear Levenberg-Marquardt least square fit we find a remarkably 
good agreement, shown in Fig. Q3] for the combined set of haloes 
of the 500/i _1 Mpc and l/i _1 Gpc simulations. It is clear that the 
three models present no substantial difference in the values of 
these distributions, meaning that the change of cosmology has no 
impact on the rotational support of the dark matter structures. 

The shape of three dimensional haloe s can be modelled 
as an ellipsoidal distribution of particles djing & Sutol 1200 2t 
lAllgood et a"i]|200fj) . characterized by the three axis a ^ b ^ c 
computed by AHF as the eigenvalues of the inertia tensor 



li.j — ^ ^ ?Ci,n3Cj,ri (23) 
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which is in turn obtained summing over all the coordinates of the 
particles belonging to the halo. 

We define the shape parameter s and the triaxiality parameter 

Tas 



5.1.4 Density profiles 

iV-body simulations have shown that dark matter hal oes can be de- 
scribe d by a Navarro Frenk White (NFW) profile dNavarro et al.l 
1 19961) . which is given by 



c 

s = — 

a 



T = 



(24) 



and we calculate the probability distributions P(T) and P(s) of 
the above parameters for all the objects above the aforementioned 
mass thresholds in our cosmological simulations, to see whether 
the VDE picture of structure formation induces changes in the 
average shape and triaxiality. Similarly to the previous case, we 
found again that halo shapes and triaxialities remain practically 
unaltered by VDE cosmolgy. This result could be expected, keep- 
ing in mind that VDE only affect background evolution: Once that 
structures start to form, detaching from the background evolution, 
they become affected by gravitational attraction only. Therefore, 
the internal structure of dark matter haloes remains generally unal- 
tered by the presence of an uninteracting form of dark energy and 
cannot be used to discriminate between alternative cosmological 
paradigms. We have also verified that these results hold also when 
taking into account different halo samples separately, i.e. the mas- 
sive ones of the l/i~ 1 Gpc simulations and the smaller belonging to 
the 500/iT 1 Mpc boxes. 



5.1.3 Unrelaxed haloes 

Before moving to the discussion of the properties of internal struc- 
ture of the haloes, and in particular the density profile, we need 
to introduce and motivate a second criterion of selection for our 
halo sample, related to the degree of relaxation of the halo. An 
additional check is necessary since only a fraction of the structures 
identified in our catalogues completely satisfies the virial condition. 
In unvirialized structures, infalling matter and merger phenomena 
may occour, heavily affecting the halo shape and thus making the 
determination of radial density profiles and concentrations unreli- 
able. In fact, unrelaxed haloes are most likely to differ from an ide- 
alized spherical or ellipsoidal shape since they have a highly asym- 
metric matter distribution, which in turn makes the det ermination of 
the ha lo ce nter an ill-defined problem, a s discussed bv lMaccio et al.l 
d2007l) and lMunoz-Cuartas et all (201 lh . Our halo finder AHF does 
not directly discriminate between virialized and unvirialized struc- 
tures giving catalogues containing both types of objects; however, it 
provides kinetic K and potential energy U for every halo identified, 
thus making the computation of the viral ra tio 2K/\U\ straigh tfor- 
ward. Following one of the criteria used bv lPrada et all d201 lh . wc 
will consider as relaxed all the haloes satisfying the condition 



2K 

W\ 



1 < 0.5 



(25) 



without introducing additional parameters. Alternative ways of 
identifying unrelaxed structures can be found throughout the lit- 
erature (e.g. iMaccio et al 1120071: iBett et al.l | 2007l: iNeto et alj|2007l; 



Knebe & Powej 120081; iPrada et all 1201 ll ; iMunoz-Cuartas et al.l 
201lMPoweretal.l2012l) : )ut since the results they give are qualita- 
tively similar for reasons of computational speed and simplicity we 
will not make use of them. The total number of haloes satisfying 
the relaxation condition is shown for every cosmology in Table [2] 



p(r) 



Po 



(26) 



where the r 3 , the so called scale radius, and the po are in principle 
two free parameters that depend on the particular halo structure. 
But po can be written as a function of the critical density as po — 
S c p c , where 



8 C 



200 



3 log(l + C )-^ 



and c = r v i T /r s is the concentration of the halo relating the virial 
radius r v (= r^oo in our case) to the scale radius r a , which will be 
discussed in detail in the following subsection. This description is 
generally valid for ACDM, but simulations of ever increased res- 
olution have actually revealed that the very central regions are not 
following the slope advocated by the NFW formula but rather fol- 
low a Sersic or Einasto profile (cf. iNavarro et al.l2004l:IStadel et al.l 
120091) . 

Here we want to check to which degree the modified cosmo- 
logical background affects the distribution of matter inside dark 
matter haloes, i.e. its density profile. All our (relaxed) objects in 
all simulations have been fitted to the Eq. (26), and to estimate the 
goodness of this fit we compute for each halo its corresponding x 2 , 
defined in the usual way 



(num) \ 2 



(num) 



(27) 



where the pi's are the numerical and theoretical overdensities in 
units of the critical density p c at the i-th radial bin and Apt is the 
numerical Poissonian error on the numerical estimate. Since differ- 
ent halo profiles will be in general described by a different number 
of radial binfl to make our comparison between different simula- 
tions and haloes consistent we need to use the reduced \ 2 



2 

Xrcd 



N, 



pts 



Ndof 



(28) 



where N pts is the total number of points used (i.e. total number of 
radial bins) and Ndof is the number of degrees of freedom (free 
parameters). 

The comparison of the distributions of the reduced \ 2 values 
for ACDM-vde, ACDM and VDE haloes belonging to the two set 
of 500/i _1 Mpc and l/i _1 Gpc simulations, shown in Fig. [T6] al- 
lows us to determine again that no substantial difference is induced 
by the VDE picture, for the same reasons discussed in the case of 
spin, shape and triaxiality distributions. The standard description 
of dark matter structures is thus not affected by the presence of a 
VDE. 



4 Note that our halo finder AHF uses logarithmically spaced radial bins 
whose number depends on the halo mass, i.e. more massive haloes will be 
covered with more bins. 
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Figure 16. Reduced \ 2 distribution for the best fit to a NFW profile, on the vertical axis we plot the total fraction of haloes whose reduced \ 2 falls within 
the horizontal axis bin value. This is shown for relaxed haloes above the 5 X 10 13 /i _1 Mq threshold belonging to the VDE-0.5, ACDM-vde-0.5 and the 
ACDM-0.5 simulations (left panel) as well as for those above the 3 X 10 14 /i -1 M Q threshold belonging to VDE-1, ACDM-vde-1 and ACDM-1 (right panel). 
The distributions show no particular difference among the three cosmologies; however, in the three lh~ 1 Gpc simulations we notice how lower resolution 
affects the \ 2 distribution, resulting in a ticker tail at higher values compared to the SOO/i^ 1 Mpc case, meaning that the fit to a NFW is on the average worse. 



5.1.5 Halo concentrations 



In the last step of the analysis of the general properties of haloes we 
will turn to concentrations, which characterize the halo inner den- 
sity compared to the outer part. This parameter is usually defined 



(29) 



where r s is the previously introduced scale radius, obtained 
through the best fit procedure of the density distribution to a NFW 
profile. We would like to remind that concentrations are corre- 
lated to the formation time of the halo, since structures that col- 
lapsed earlier tend to have a more compact center due to the fact 
that it has more time to accrete matter from the outer parts. Dy- 
namical dark energy cosmologies generically imply larger c val- 
ues as a consequen ce of earlier structu r e formation, as found in 
wor ks like those by iDolag et al. I d2004l) . I Bartelmann et al. 
and lGrossi & Springell ( 20091) . In fact, since the presence of early 
dark energy usually suppresses structure growth, in order to repro- 
duce current observations we need to trigger an earlier start of the 
formation process, which on average yields a higher value for the 
halo concentrations. However, this result does not hold in the case 
of coupled dark energy, where the increased clustering strength in- 
duced by a fifth force sets a l ater start of structure formation, as 
discussed in iBaldi et alj JIoToh . 

In the hierarchical picture of structure formation, concentra- 
tions are usually inversely correlated to the h alo mass as more 
massive objects fo r m later; A/-body simu l ations dDolag et al J 20041 ; 



Prada 
dCorr 



Prada et al.l l201ll; Munoz-Cuartas e t al.l 1201 lh and observations 



Comerford & Nataraianll2007l;IOkabe et al.ll2010l;lsereno & Zitrinl 



120121) have in fact shown that the relation between the two quanti- 
ties can be written as a power law of the form 



logc = a(z)log( 



M-2 



00 







+ b(z) 



(30) 



where a(z) and b(z) can have explicit pa rametrizations as 
functions of redshift and cosmology (see e.g. iNeto et all 120071 ; 
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Figure 17. Best fit of the mass-concentration relation for the combined sam- 
ple of relaxed haloes belonging to all the ACDM, ACDM-vde and VDE 
simulations. The points represent the average concentration values for the 
relaxed haloes in the corresponding mass bin; circles are for ACDM, trian- 
gles for ACDM-vde and squares for VDE. Empty dots stand for bins de- 
termined using haloes belonging to the 1 h~ x Gpc simulations; filled ones 
refer to the 500 Mpc ones. The Poissonian error bars are computed 
using the number of selected haloes within each mass bin. 



IPrada et al.ll2oTT| ; iMufioz-Cuartas et alj|201 ll) . We can use our se- 
lected halo samples at z = from the 500 h" 1 Mpc and 1 fe~ 1 Gpc 
simulations to obtain the a(z — 0) and b(z — 0) values for the 
ACDM, ACDM-vde and VDE cosmologies; the results of the best 
fit procedure to Eq. d30t are shown in Table [3] 

These values are in goo d agreement with the ones f ound 
by, for instance. IDolag et al.l d2004l) . iMaccio et all J2008h and 
IMufioz-Cuartas et al.l J201 lh (who quote for ACDM values of 
a{z = 0) « -0.097 and b « 2.01); the ~ 10% discrepancy 
observed with their results is due to the fact that our results are 
obtained over a smaller mass range, 5 x 10 13 -2 x 10 15 /i _1 M , 
whereas the previously cited works study it over an interval larger 
by more than three orders of magnitude, 1O 1o -1O 15 /i -1 M . Still, 
according to our results, the c — Al relation for both the VDE and 
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Table 3. Best fit values for the mass-concentration relation for z = 0, ob- 
tained fitting the relation Eq. < |30j to the relaxed haloes concentrations and 
the two point correlation function to a power law (ro/r) 7 for the ACDM, 
ACDM-vde and VDE cosmologies, ro values are given in h _1 Mpc. 



Model 


o 


b 


'■(i 


7 


ACDM 


-0.115 


2.11 


13.4 


-1.79 


ACDM-vde 


-0.112 


2.21 


12.1 


-1.91 


VDE 


-0.098 


2.17 


10.1 


-1.94 



The identification of "sister haloes" among the different cos- 
mologies can be done using the AHF tool MergerTree, which 
determines correlated structures by matching individual particles 
IDs in different simulation snapshots. For a more elaborate discus- 
sion of its mode of op eration we refer the reader to Section 2.4 in 
Libesk ind et al. where it has been described in greater de- 

tail. This time we decided to restrict our halo sample further only 
picking the first 1000 most massive (ACDM) haloes. The criterion 
of halo relaxation has of course also been taken into account when 
dealing with profiles and concentrations. 




0.01 r 



0.001 



Figure 18. Distribution of the ratio between the actual concentration and 
the expected one (cf. Eq. )30t ) and its fit to a lognormal distribution. 



ACDM-vde case is characterized by a shallower a exponent and 
a larger b. Although the magnitude of these changes is different 
in the two models, we can safely conclude that also in this case 
the results are mainly parameter-driven, i.e. due to the larger value 
of Qm- Furthermore, the large error bars for M > 10 15 7i _1 M 
scales, due to the low statistics of massive haloes complying the 
relaxation requirements, makes it difficult to determine to which 
extent the differences in the best fit relations among ACDM-vde 
and VDE could be eventually reduced in the presence of a larger 
sample. 

We also need to mention that in our simulations the actual halo 
concentrations do not precisely follow equation Eq. ( 130) but rather 
scatter around it, as can be seen in Fig. [T7] where the average c 
per mass bin is plotted against the corresponding best fit relations. 
This is not really surprising, s ince observations (ISereno & Zitrinl 
120121) and TV-body simulations l lDolag et al.l2004l) have shown that 
halo concentrations are lognormally distributed around their theo- 
retical value calculated using Eq. i30) . In Fig.[T8]we show that this 
is indeed the case: the distribution of the c(M)/cfu(M), where 
Cfu(M) is the theoretical concentration value for a halo of mass 
M, extremely close to a lognormal one with an almost model inde- 
pendent dispersion a ~ 0.4. 



5.2 Cross Correlation 

The next step in our analysis consists of studying the properties of 
the (most massive) cross correlated objects found in the three mod- 
els at z = 0. Whereas in the previous section our focus was on the 
distribution of halo properties, this time we aim at understanding 
how they change switching from one model to another. 



5.2.1 Mass and spin parameter 

In the two upper panels of Fig.[l9]we show the ratios of the masses 
M and spin parameter A for all the cross correlated sets of simula- 
tions; in each panel we show the ratios for the 500/i _1 Mpc simu- 
lation boxes while on the right the l/i~ 1 Gpc ones. Both VDE and 
ACDM-vde show average mass and spin values scattered around 
values larger than one when compared to ACDM, whereas the cross 
comparison of VDE to ACDM-vde shows average ratios close to 
unity at all mass scales. This substantial increase in the ratios is 
due to the earlier beginning of structure formation, triggered by the 
larger Q.m and 0$, as the comparison VDE/ ACDM-vde shows. As 
we already did in Section |5~T] when looking at the halo properties 
in general, we conclude that also when observing the same halo 
evolved under different cosmologies, the main effects are deter- 
mined exclusively by the set of cosmological parameters chosen, 
being the imprint of the cosmological background evolution sub- 
stantially negligible in this case. This makes the identification of a 
cosmic vector through the determination of halo properties impos- 
sible, since the background dynamics, which distinguishes VDE 
from any other non interacting dynamical dark energy model, does 
not leave any observable imprint on these scales. 



5.2.2 Halo concentrations and internal structure 

As done in the previous section, in the determination of halo pro- 
files and concentrations properties we discard unrelaxed haloes, but 
this time in a way so that our halo sample will still be composed of 
the first 1000 haloes satisfying condition Eq. ( 125) . This same halo 
sample has been used also in the study of the M V i r — Zf orm , in 
order to be able to compare these results with the ones obtained 
from concentrations consistently - although in principle formation 
redshifts are well defined even for unrelaxed haloes. Again, our 
procedure consists in fitting all the selected structure to a NFW 
profile, from which we will be able to derive the concentration pa- 
rameter c and a measure for the quality of the fit x 2 ; we will then 
compare these results in each cross identified objects to see how a 
given halo structure changes when evolved under a different cos- 
mology. Although not shown here, no particular trend in the dif- 
ferences among ACDM, ACDM-vde and VDE pictures has been 
be found for either NFW yj 2 , shape and triaxiality; since in all the 
cases the ratios of these properties among cross correlated haloes 
are centered around unity. Not surprisingly, we also find again a 
generally higher average value for the concentrations in VDE and 
ACDM-vde with respect to ACDM, (see Fig. [19} a result which 
again can be explained by the larger value of Qm and erg. Similar 
concentrations for VDE and ACDM-vde haloes, shown in the up- 
per right panel of Fig.[l9]can be also understood as a consequences 
of the similar masses of the haloes examined and the similar c — M 
relations found for the two cosmologies. However, even if from 
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Figure 19. Mass and spin parameter correlation ratios for the first 1000 (relaxed and possibly unrelaxed) haloes. Panels on the left show the VDE/ACDM 
results, the ones in the center to ACDM-vde/ACDMwhile on the right the ratio VDE/ACDM-vde is plotted. Cross identified objects are characterized by larger 
masses in VDE and ACDM-vde as a consequence of the higher Qj^ and <rg normalization value. 



ACDM-vde cosmology we conclude that the different choice of 
f!ju, can explain in this case higher halo concentration, we need to 
remind that such a results is also a general feature of the dynami- 



cal na t ure of the dark energy flu id, as already found by Dolag et all 
foil 4 . 1 Bartelmann et al. | [2006h and lGrossi & Springelld2009h . 



6 CONCLUSIONS 

In this work we presented an in-depth analysis of the results of a 
series of A^-body dark matter o nly simulations of the Vector Dar k 
Energy cosmology proposed bv lBeltran Jimenez & Marotol ( l2008l) . 
The main emphasis has been on the comparison to the standard 
ACDM paradigm, using a mirror simulation with identical number 
of particles, random seed for the initial conditions, box size and 
starting redshift. An additional series of simulations for a ACDM- 
vde cosmology have also been run using the VDE values for £Im 
and as within a standard ACDM picture, to disentangle the effects 
of the parameter induced modifications to the dynamical ones com- 
ing directly from the VDE model. 

The use of a modified version of the GADGET-2 code required 
us to check the results with particular care. A consistency check of 
our simulations was performed by comparing the numerical results 
for the evolution of the growth factor to the analytical calculations, 
finding an excellent agreement between the two. We further had to 
adapt the halo finding procedure, due to the fact that the critical 



density as a function of redshift p c (z), entering the definition of 
the halo edges, takes different values in VDE. Once halo catalogues 
had been obtained, we carried our analysis at two different levels, 
namely: 

• we studied the very large scale clustering pattern through the 
computation of matter power spectra, mass, void, and two-point 
correlation functions; 

• we analyzed halo structure, comparing statistical distribu- 
tions and averages of spin parameters, concentrations, masses and 
shapes. 

In the first point, making use of the full set of simulations, 
our analysis covered the whole masse range 10 12 -10 15 /i^Mq 
as well as different redshifts, so that we could make specific VDE 
model predictions for the number density evolution n(> M, z) and 
growth index 7(2). A distinctive behaviour, very far from the stan- 
dard ACDM results, has been found for 7(2), and, in particular, for 
the mass function that in VDE cosmology can be up to 10 times 
larger than the standard ACDM one. The latter result is due to the 
earlier onset of structure formation and we have mentioned how it 
can be used to address current ACDM observational tensions with 
large cl usters at z > 1 and possibly with early reionization epoch 
(cf. also lCarlesi et alj|20lj) . At z = 0, we have shown that VDE 
cosmology requires the standard set of parameters entering the Tin- 
ker mass function to be modified; however, since the investigation 
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of the linear perturbation theory is still ongoing, we lack the instru- 
ments to shed more light on this aspect. 

Computing the cumulative mass function at different redshifts 
and making use of the ACDM-vde simulations we have also ob- 
served how the condition Hvde(z) > Hacdm(z), holding up to 
z » 1, induces a relative suppression of structure growth in this 
cosmological model, an effect that clashes with the increased mat- 
ter density and erg. In fact, while on the one hand higher values of 
these parameters enhance the formation of a larger number of ob- 
jects, on the other hand, background dynamics suppresses cluster- 
ing and growth. The interplay and relative size of these effects has 
been studied using the ACDM-vde simulations, showing that, for 
example, faster expansion in the past determines for VDE an expec- 
tation of clusters with M > 10 14 /i _1 Mq up to « 5 times smaller 
than what a simple increase in as and Qm would determine. This 
effect has been also seen in the void distribution, where suppression 
of clustering prevents small structures to merge into larger one and 
to rather spread in the field, so that underdense regions happen to 
be smaller and rarer than in ACDM and ACDM-vde. In these latter 
cosmologies, in fact, a higher contrast between populated and less 
populated regions is observed both in the power spectrum and in 
the colour coded matter density. 

In the second part of our work we have focused on the study 
of internal halo structure. We found that VDE cosmology does 
not induce deviations in the functional form of the dark matter 
halo density profile s, which are still well described by a NFW 
dNavarro et al.ll 1996b profile, nor in the distributions for the con- 
centrations and spin parameters, which are of the lognormal type 
as in ACDM. Shape and triaxiality are also unaffected: the distri- 
butions for the relative parameters are identical and peaked in at the 
same values in all the three cosmologies. The above results are a di- 
rect consequence of the fact that dark matter haloes, once detached 
from the general background evolution driven by the cosmic vec- 
tor, evolve by means of gravitational attraction only; which is unaf- 
fected by the specific nature of dark energy. A net effect can be seen 
in masses, whose average values tend to be on the larger than in the 
ACDM case by a factor of ~ 2, a straightforward consequence of 
the larger Q,m and as, as can be shown by a direct comparison of 
VDE to ACDM-vde results, that turn out extremely close in these 
cases. On the other hand, the different background evolution seems 
to affect c — M relations only slightly, changing the power law in- 
dex a(z) and normalization b(z) by a 15%. In this case we have 
also found that these values in general agree with previous result s 
from early dark energy studies such as those bv lDolag etaTl(l2004l) , 
even though in this case it would certainly be necessary to test the 
relation down to smaller mass scales, where a better tuning of the 
parameter would be also possible, and with a larger statistics on 
the higher scales. However, in general, most of the halo-level ef- 
fects which seem to characterize VDE can be simply explained in 
terms of the different cosmological parameters, as we did compar- 
ing these results to the outcomes of ACDM-vde simulations. For 
the first time then, through the results of the series of A^-body sim- 
ulations, we have shown that VDE cosmology provides a viable 
environment for structure formation, also alleviating some obser- 
vational tensions emerging with ACDM. We have seen how the 
peculiar dynamics of this model leaves its imprint on structure for- 
mation and growth, and in particular, how it affects predictions for 
large scale clustering and halo properties. However, a close com- 
parison of the deep non-linear regime results with different sets of 
observational data still needs to be performed, challenging us to im- 
prove the accuracy of our simulations and at the same time devise 
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new and reliable tests which may shed some light not only on VDE 
but on the nature of dark energy in general. 
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